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ASSESSMENT  OF  DAMAGE  AGGREGATION  FOR  WEAPONS  SALVOS 


Toke  Jayachandran  and  I.Bert  Russak 


1 .  INTRODUCTION.  Estimation  of  the  expected  proportion  of  an  area  target 
damaged  by  a  salvo  of  weapons  is  an  important  problem  in  weapons 
systems  effectiveness  studies.  Analysts  sometimes  use  an  approximate 
extant  formula,  referred  to  herein  as  the  “empirical  rule”,  for  this 
purpose.  The  genesis  and  the  rationale  underlying  this  formula  are  not 
readily  available.  Esary  [  1  ]  ,  [  2  ]  introduced  a  “proportional  effects 
model”  for  which  the  empirical  formula  is  a  first-order  approximation  to 
the  expected  cumulative  damage  to  the  target.  Assessment  of  this 
damage  exactly  would  require  a  knowledge  of  the  probability  distribution 
of  the  number  of  weapons  from  the  salvo  that  will  impact  on  the  target. 
Esary  investigated  the  accuracy  of  the  empirical  rule  approximation  for 
three  choices  of  the  weapons  impact  distribution.  His  conclusion  was  that 
the  empirical  rule  consistently  overestimates  the  cumulative  damage  for 
all  possible  impact  distributions  and  that  the  approximation  can  be  quite 
poor  for  many  realistic  situations.  The  advantage  in  using  the  empirical 
rule  is  that  one  can  compute,  approximately,  the  cumulative  damage  using 
only  the  mean  or  expected  number  of  weapons  that  will  hit  the  target, 
without  knowing  the  actual  distribution  of  hits.  In  this  report,  we 
propose  modifications  to  the  empirical  rule  that  use  in  addition  to  the 
mean,  the  variance  of  the  hit  distribution.  These  modifications  result  in 
significantly  improved  measures  of  the  expected  cumulative  damage  to 
the  target. 

2.  PROBLEM  DFRCRIPTIQN.  A  salvo  of  n  weapons  is  launched  at  an  area 
target.  The  number  of  weapons  that  will  actually  hit  the  target  is  a 
random  variable  N  whose  possible  values  are  0,1,2,.  .  .,n.  Each  hit  will 
result  in  damaging  a  certain  proportion  of  the  target.  The  cumulative 
damage  to  the  target  is  then  a  random  variable  D  which  depends  on  N  the 
random  number  of  weapons  that  impact  the  target.  The  quantity  of 
interest  is  E(D)  the  expected  (or  mean)  cumulative  damage.  To  estimate 
E(D)  exactly,  the  probability  distribution  of  N  the  number  of  weapons  that 
impact  the  target,  the  target  configuration  and  the  damage  aggregation 
scenario  for  multiple  hits  are  needed.  However,  in  practice  analysts  are 
known  to  use  the  empirical  rule 
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E(D)  =  1  -  (1-cl)E(N) 


(1) 


where  d  is  the  damage  proportion  for  a  single  hit  and  E(N)  is  the  expected 
number  of  weapons  out  a  salvo  of  n  weapons  that  hit  the  target.  As 
mentioned  earlier,  the  genesis  of  this  formula  is  not  well  documented. 
Esary  [  1  ]  considered  the  following  “proportional  effects”  model  that  can 
be  used  to  provide  a  rationale  for  the  empirical  rule.  The  premise  of  this 
model  is  that  a  proportion  d  of  a  pristine  target  is  damaged  by  the  first 
hit  and  each  additional  hit  damages  the  same  proportion  d  of  that  part  of 
the  target  not  previously  damaged.  Then,  it  can  be  shown  that  if  N 
weapons  hit  the  target,  the  cumulative  damage  is 

D(N)  =  1  -  {1-d)N  (2) 

and  E[D(N)]  ,  or  simply  E(D),  is  the  expected  cumulative  damage  to  the 
target.  Replacing  N  with  E(N)  in  D(N)  will  result  in  the  empirical  rule  E(D) 
which  will  be  different  from  E(D),  in  general.  To  evaluate  E(D)  exactly, 
the  probability  distribution  of  N,  the  number  of  weapon  hits  is  needed. 
Esary  proved  that  the  empirical  rule  consistently  over  estimates  the 
expected  cumulative  damage  i.e.,  E(D)  >  E(D)  for  every  possible  hit 
distribution.  If  it  is  assumed  that  the  distribution  is  Binomial  i.e.,  if  each 
of  the  n  weapons  hits  the  target  independently  with  the  same  probability 
p,  the  expected  cumulative  damage 

Eb(D)  =  1  -  (1-pd)n  (3) 

and  the  estimate  using  the  empirical  rule  for  this  case  is 

Eb(D)  =  1-(1-D)np  (4) 

For  a  Poisson  (sometimes  used  as  an  approximation  to  the  Binomial)  hit 
distribution  with  parameter  X=np, 

Ep(D)  =  1  -  e-^d  (5) 

and  Ep(D)  will  be  the  same  as  for  the  binomial  distribution.  The 
corresponding  formulas  for  a  Uniform  hit  distribution  are 

Eu(D)  =  1  -  (1-d)n+i/(n+1)d  (6) 
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Eu(D)  =  1  -  (1-d)n/2. 


(7) 


For  the  Binomial  and  the  Poisson  distributions,  the  exact  values  of  E(D) 
were  compared  with  the  corresponding  approximate  values  for  the 
empirical  rule,  for  various  choices  of  n,  p  and  d,  to  demonstrate  the 
magnitudes  of  the  errors  in  the  approximation.  The  resulting  graphs  are 
presented  in  figures  1-16.  It  can  be  seen  that  the  two  measures  are 
reasonably  close  for  small  values  of  d  (0<d<.5)  and  they  diverge  markedly 
as  d  approaches  1  for  all  choices  of  n  and  p.  However,  the  disparity 
decreases  as  n  the  number  of  weapons  in  the  salvo  increases. 

3.  A  BECOND-ORDER  TAYLOR  APPROXIMATION.  The  Taylor  approximating 
polynomial  of  degree  two  about  the  point  X=|j.x  for  a  function  f(X)  of  a 
random  variable  X  with  |j.x=  E(X)  is 

f(X)  =  D(  Hx)  +  f'(Hx)  (X-  Hx)  +  f■(^*)(X-^x  )^/2  (8) 

and  taking  expectations  on  both  sides 

E[f(X)]  =  f(nx)  +  %x)o^x/2  (9) 

where  o2x  =  V(X)  the  variance  of  the  random  variable  X.  Applying  this 
result  to  D(N)  in  equation  (2)  above,  we  get  the  second-order  Taylor 
approximation  to  E[D(N)], 

T[D(N)]  =  D(pn)  +  D”(Pn)  /2  (10) 

=1-{l+[ln(1-d)]2a2N /2}(1-d)^lN 

where  In  is  the  natural  logarithm  function  and  Pn  ,  n  are  the  mean  and 
the  variance  of  the  random  variable  N.  We  propose  T[D(N)],  the  Taylor 
approximation  as  a  better  alternative  to  the  empirical  rule  In  equation 
(1).  Whereas  only  pn  is  needed  to  use  the  empirical  rule,  both  pn  and 
are  required  for  the  use  of  this  second-order  approximation. 

Graphs  of  the  Taylor  approximation  are  included  in  figures  1-16;  they 
show  significant  improvements  in  the  approximation  relative  to  the 
empirical  rule,  for  all  choices  of  the  parameters  n,p,d.  However,  for 
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larger  values  of  d  {d>.8)  the  Taylor  approximation  errors  are  still  a  little 
on  the  high  side  (albeit  smaller  than  that  for  the  empirical  rule), 
especially  for  small  values  of  n  the  number  of  weapons  in  the  salvo.  To 
remedy  this  problem,  we  propose  an  ad  hoc  correction,  determined  through 
trial  and  error,  to  the  Taylor  approximation,  to  further  reduce  the  error  at 
the  high  values  of  d.  Specifically,  we  suggest  the  modified  formula 

M[D(N)]  =  T[D(N)]  -  D4(1-P2)/2(N+1)  .  (11) 

Approximation  error  plots  for  the  modified  formula  are  also  shown  in 
figures  1-16.  The  modification  appears  to  work  well,  in  most  cases,  as 
can  be  seen  from  the  closeness  of  the  graphs  to  those  for  the  exact 
formula.  The  exceptions  are  for  the  case  of  the  Uniform  hit  distribution 
which  is  a  flat  distribution  i.e.,  its  graph  is  rectangular  in  shape.  The 
Uniform  distribution  as  a  model  for  the  number  of  hits  in  a  salvo  of  n 
weapons  is  not  very  realistic  since  it  would  imply  that  the  probability  of 
0  hits  or  1  hit  or  2  hits  ...  or  n  hits  is  exactly  the  same.  The  distribution 
was  included  in  this  study  to  indicate  the  behavior  of  the  approximations 
in  the  case  of  a  relatively  flat  hit  distribution. 

Error  percentages  (100x(approximate  -  exact)/exact  )  are  presented  in 
Tables  1-3  for  each  of  the  three  approximation  formulas  viz..  Empirical, 
Taylor  and  Modified  for  n=4,6,8;  p=.4,.6,.8:  and  the  three  hit  distributions, 

Binomial,  Poisson  and  the  Uniform.  These  tables  show  more  clearly  the 
improvements  attained  by  using  the  approximations  proposed  in  this  paper. 


4.  CONCLUSIONS.  As  shown  by  Esary,  the  empirical  rule  grossly 
overestimates  the  expected  cumulative  damage.  We  proposed  two  new 
approximate  rules  viz.,  the  Taylor  approximation  and  the  Modified 
approximation  that  consistently  provide  better  results  than  those  for  the 
empirical  rule.  The  modified  approximation  provides  the  best  measures 
unless  the  hit  distribution  is  somewhat  flat,  i.e,  the  distribution  is  close 
to  the  uniform  distribution;  in  which  case  the  Taylor  approximation 
appears  to  be  the  best  choice. 
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TABLE  I 

PERCENT  ERRORS  FOR  THE  BINOMIAL  DISTRIBUTION 
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TABLE  II 

PERCENT  ERRORS  FOR  THE  POISSON  DISTRIBUTION 
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TABLE  III 

PERCENT  ERRORS  FOR  THE  UNIFORM  DISTRIBUTION 
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